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I. Introduction 

Collective phenomena in systems far from thermal equilibrium have been a subject of 
extensive studies in recent years. Usually these systems are driven out of equilibrium by 
external fields, such as electric field in the case of conductors, pressure gradient in the case 
of fluids, temperature gradient in the case of heat conductors, chemical potential gradient in 
the case of growth problems and many others . These driving fields are very common in 
nature and are found in a large variety of physical systems such as granular and traffic flow 
1^-^, gel electrophoresis @J^, superionic conductors [II^jO] to give a few examples. In many 
cases these systems reach a steady state, which unlike the equilibrium case, is characterised 
by non-vanishing currents. In these lectures we consider possible collective phenomena and 
phase transitions which may take place in such steady states. The main problem in studying 
nonequilibrium systems is the lack of general theoretical framework within which they could 
be analysed. As a result they are far less understood as compared with equilibrium systems 
where the Gibbs picture provides such a theoretical framework. 

Before discussing nonequilibrium systems it is useful to consider briefly systems in ther- 
mal equilibrium. Here decades of studies have yielded a fairly detailed understanding of 
their thermodynamic behaviour. Many rules which govern phase transitions occurring in 
these systems have been derived. For example, it has been shown that the critical exponents 
associated with a phase transition may be classified into universality classes. These classes 
do not depend on the detailed interactions in the system but rather on a few parameters 
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such as the symmetry of the system and of the order parameter associated with the tran- 
sition, the dimensionahty of the system and the range of interactions. Therefore, in order 
to study theoretically the critical behaviour of a given system it is sufficient to analyse the 
simplest possible model which belongs to the same universality class. For reviews see, for 
example, |T2|JT^. It has also been shown that phase transitions and spontaneous symmetry 
breaking do not take place at low dimension. In particular, no phase transition is expected 
to take place in a one dimensional system at finite temperatures as long as the interactions 
are short range Moreover, breaking of continuous symmetry may take place under the 



same conditions only in dimensions higher than two Other rules derived by Landau 

relate the nature of the transition, namely whether it is first order or continuous, to the 



symmetry of the systems . If the symmetry allows a third order term in the expansion of 



the free energy in the order parameter, such as in the case of the transition from a liquid to a 



nematic liquid crystal phase [|T6[, the transition cannot be continuous and is necessarily first 
order. On the other hand if the symmetry is such that no third order term is allowed, such 
as in the transition from a paramagnetic to a ferromagnetic phase, the transition may either 
be first order or continuous, depending on the details of the interactions. The Gibbs phase 
rule is another very useful example of a rule which governs the phase digrams of systems in 
equilibrium [|l4l . It deals with fiuids composed of c components. The thermodynamic phase 



space of such systems is of c + 1 dimensions, associated with the temperature, pressure, and 
c — 1 chemical potentials. According to the rule, the manifold in this space on which n 
different phases coexist isofD = 2 + c — n dimensions. Another rule deals with the phase 
diagram near a triple point, where three coexistence lined meet. According to this 180° rule, 
each of the three angles defined by the intersecting coexistence lines must be less than 180°. 



This is a direct result of the convexity of the free energy |T7 



These rules and many others, some of which are related to disordered systems ||T8|,|T9 
provide extremely useful tools for analysing and understanding phase diagrams and critical 
behaviour of models and physical systems in equilibrium. By simply identifying the symme- 
try of the system and the nature of the order paramenter involved in the phase transition 



one can usually find the universality class of the transition and even obtain a rough idea of 
the possible phase diagram. 

Our degree of understanding of collective behaviour far from thermal equilibrium is at a 
much more primitive stage. Since a general theoretical framework for studying nonequilib- 
rium phenomena does not exist, one cannot derive similar rules which would be as general 
as those for equilibrium systems. Rather, one has to resort to studying specific models and 
probe the resulting types of phase diagrams and phase transitions, with the hope that some 
general picture might emerge. 

In the present lectures we consider stochastic driven systems in one dimension and discuss 
some interesting collective behaviour which they display. Unlike equilibrium one dimensional 
systems which do not exhibit phase transitions, non equilibrium systems exhibit a rich variety 
of collective phenomena such as first order and continuous phase transitions, spontaneous 
symmetry breaking (SSB), phase separation, slow coarsening processes and many others. 
Mechanisms which lead to these phenomena are discussed. 

The article is organised as follows: in Section |I| the concept of detailed balance is dis- 
cussed, the lack of which is characteristic of nonequilibrium systems. A necessary and 
sufficient condition for the existence of detailed balance is presented. In Section |TT| a simple 
driven model, the totally asymmetric exclusion process, is introduced and its phase diagram 
for a system with open boundaries is calculated using a mean field approximation. The 
phase diagram exhibits several phases separated by first order and continuous transitions. 
The matrix method which enables one to obtain exact results for steady state properties 



is outlined in Section A model which displays spontaneous symmetry breaking in one 
dimension is introduced in Section ^ and a model exhibiting phase separation accompanied 
by slow coarsening processes is described in Section [V^. Open problems and perspectives 
are briefiy discussed in Section |V11. 
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II. Detailed balance and driven systems 



In this section we make some general considerations concerning the evolution of dynam- 
ical systems. Let C be a microscopic configuration, and let P{C, t) be the probability that 
the system is in the microscopic configuration C at time t. The dynamics of the system 
is defined in terms of the transition rates W{C C) from a configuration C to C . The 
equation which governs the evolution of the distribution function P(C, t) takes the form 



The first sum represents the rate of fiow, in configuration space, of probability into C while 
the second sum corresponds to the outgoing fiow from this configuration. In a steady state 
the two terms are equal, yielding zero net fiow from any configuration. 

Systems in thermal equilibrium are characterised by an energy function, or a Hamilto- 
nian, E{C). The steady state distribution P{C) is proportional to e~^^^^^''^'^ , where T is 
the temperature and ks is the Boltzmann constant. Given an energy function E{C) one can 
always find transition rates W{C C), such as the Metropolis rates, which obey detailed 
balance. Here the two sums cancel term by term 



for any pair of configurations C and C. 

On the other hand dynamical systems are not defined by an energy function but rather 
by transition rates. When a system is not in thermal equilibrium, the resulting steady is 
such that detailed balance @ is not satisfied. We will basically use this lack of detailed 
balance definition of nonequilibrium. 

Given the dynamics of a system, namely the transition rates, it is of interest to know 
whether or not detailed balance is satisfied. Since, in general, the steady state distribution 
cannot be calculated, a direct check of the detailed balance condition (0) is not possible. 
Thus a criterion for existence of detailed balance which is based directly on the transition 
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di 



w[c' c)p{c', t)-J2 ^ic c')p{c, t). 



c c 



(1) 




(2) 



4 



rates and does not require the knowledge of the steady state is highly desirable. Such a 
criterion is provided by the following equations. Let Ci, C2, . . . , Cfc be a set of microscopic 
configurations. A necessary and sufficient condition for the existence of detailed balance is 
that for any such set one has 

W{1 2)W{2 ^ 3) . . . W{k ^ 1) = W{1 k)W{k ^ k ~ 1) . . . W{2 1), (3) 

where for simplicity we have denoted Ci hj i. It is easy to check that (H) is a necessary 
condition. When detailed balance is satisfied one may replace W{i ^ i + 1)/W{i + 1 — > i) 
by P{i + 1)/P{i), where P(i) is the steady state distribution with respect to which detailed 
balance is satisfied. Using these relations @ is easily verified. 

To demonstrate that this is a sufficient condition as well, we use (|^) to derive the steady 
state distribution. We start with an arbitrary configuration 1 and denote its steady state 
weight by -P(l). The weight of states 2 which are directly connected with 1 (namely, for 
which W{1 —>■ 2) > 0), may thus be defined using the detailed balance relation, P{2) = 
P{1)W{1 2)/W(2 1). This process may then be repeated to define the weights of 
states directly connected with states 2 etc, until all microscopic configurations have been 
reached. The weight of a microscopic configuration k which may be reached from 1 via 
intermediate states 2, 3, . . . , /c — 1 is thus given by 

W(1^2)...W(k-l^k) 

'^'•^>w{k^k-i)...w(2^iy ^' 

For this procedure to be self-consistent one has to verify that any path between configurations 
1 and k yields the same P{k). It is a straightforward matter to show that this follows directly 
from d^). 

Thus, to demonstrate that a dynamical system defined by its transition rates is not in 
thermal equilibrium it is sufficient to find a single path in configuration space for which 
(0) is not satisfied. This is usually quite easy to check, making it a very useful criterion. 
When is not satisfied the system exhibits non-vanishing probability currents between 
configurations which is an indication of the system not being in thermal equilibrium. 



A simple prototypical model of driven systems, termed the 'standard model', was intro- 
duced by Katz et al pO|,p^. This is a driven lattice gas model defined on a hypercubic lattice 
with periodic boundary conditions. Each site i is either occupied by a particle or is vacant, 
with cxj = 0, 1 being the occupation number. In the absence of drive, an Ising Hamiltonian 
is assumed 

= -J^o-idj, (5) 

{ij) 

where the sum is over nearest neighbour (nn) sites {ij). The evolution of the system is defined 
by Kawasaki dynamics, allowing for particles to hop between nearest neighbour sites. Let 
C and C be two configurations obtained from each other by an interchange of a single pair 
of nn occupation numbers (Ti and o"-,-. The transition rate between C and C may be taken 
as the Metropolis rate W{C C) = w{(3AH), where AH = H{C') - H{C) and w{x) = 
min(l,e~^). This dynamics leads to the expected Boltzmann equilibrium distribution. In 
d > 1 dimensions the system exhibits the usual Ising transition from a homogeneous phase 
at high temperatures to a phase separated state at low temperatures. 

Introducing a driving field E along one of the axes, the transition rates are modified by 
adding a term uE to AH where 

u = -1 , , +1, (6) 

for a hop along, transverse or opposite to the field direction, respectively. The transition 
rates are thus given by 

W{C ^ C) =w{/3{AH + uE)). (7) 

Due to the periodic boundary conditions in the direction of the driving field, these rates do 
not obey detailed balance, and the steady exhibits non-vanishing currents. 

In spite of the simplicity of the model, no exact results for the steady state properties 
are available (exact in one dimension). Extensive numerical studies of this model in two 
dimensions demonstrate that the phase separation transition which exists in zero drive, 
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persists for non-zero driving fields. At low temperatures the system exhibits stripes of high 
density and low density regions which are oriented along the field direction. These stripes 
coarsen with time leading to a phase separated state. 



III. Asymmetric exclusion process in one dimension 

In this section we consider the phase diagram of the one dimensional 'standard model' 
defined in the previous section for J = 0. Here the only interaction between the particles is 
the hard-core interaction which prevents more than one particle from occupying the same 
site. This process is called asymmetric simple exclusion process (ASEP). Furthermore, we 
consider the limit E oo, called totally asymmetric simple exclusion process (TASEP). In 
this limit particles are restricted to move only to the right, with no backward moves. This 
model turns out to be sufficiently simple to allow for exact calculation of some of its steady 
state properties. In spite of its simplicity, the model with open boundary conditions exhibits 
a rather rich and complicated phase diagram, displaying both continuous and discontinuous 
phase transitions (see below). This is clearly a direct consequence of the fact that the 
dynamics is a nonequilibrium one. The model and many variants of it have been a subject 



of extensive studies in recent years [pTH30 



The dynamics of the model is defined as follows: at any given time a pair of nn sites is 
chosen at random. If the occupation numbers of these sites are (+ 0) an exchange is carried 
out 

+ 0^0+, (8) 

with rate 1. All other configurations, remain unchanged. Here and in the following we 
interchangably use 1 or + to denote an occupied site. For periodic boundary conditions 
the system reaches a trivial steady state in which all microscopic configurations have the 
same weight. This may be verified by direct inspection of the master equation dTj). It is 
easy to see that the number of configurations to which a given configuration C may flow 



7 



is equal to the number of configurations flowing into C. To verify that this is the case 
note that a microscopic configuration C is composed of alternating segments of +'s and O's. 
According to the dynamics (§) C may be exited when the rightmost particle in one of the + 
segments moves one step to the right. Thus the number of configurations C may flow into 
is equal to /, the number of + segments in this conflguration. Similarly C may be reached 
when the leftmost particle in one of the + segments hops into its position. The number of 
configurationd fiowing into C is therefore also equal to /. Since all non-vanishing transition 
rates are 1, the state where all configurations have equal weights is stationary. Therefore 
in the steady state the system exhibits no correlations, apart from the trivial correlations 
arising from the fact that the overall density of particles is fixed. 
The steady state current J is given by 

J=(a,(l-a,+i)), (9) 

where the brackets denote a statistical average with respect to the steady state weights of 
the microscopic configuration. Since in the steady state the system exhibits no correlations 
the current may be written, in the large system limit, as 

J = p{l-p), (10) 

where pi = (cXi) is the density at site i, and the index i is omitted in ( p!o|) since the average 
density is homogeneous, independent of i. Equation (p!0D , relating the current to the density 
is known as the fundamental relation (or fundamental diagram). The interesting feature 
in this relation is that the current is not a monotonic function of the density but rather it 
exhibits a maximum at p = 0.5. This feature is a result of the hard-core interaction between 
the particles and it affects rather drastically the steady state properties of the system when 
open boundary conditions are considered. 

We now turn to the model with open boundary conditions. Here, particles are introduced 
into the system at the left end, they move through the bulk according to the conserving 
dynamics (^), and leave the system at the right end. To be more specific, at the left 
boundary {i = 1) the move 
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FIG. 1. A schematic representation of the totally asymmetric exclusion process with open 
boundary conditions. 

0^+, (11) 

is carried out with a rate a. Similarly, at the right boundary [i = N) one takes 

+ - 0, (12) 

with a rate f3. For a schematic representation of the model see Figure |^. The overall 
dynamics is non-conserving. Particles are conserved in the bulk but are not conserved 
at the boundaries. Unlike the case with periodic boundary conditions, the steady state 
distribution is not trivial, and correlations between the densities at different sites do not 
vanish. However, far from the boundaries the distribution function is expected to be well 
approximated by the homogeneous one, suggesting that correlations are small. 

We are interested in the steady state of this model for given rates a and (3. For large 
a and small (3, namely for a large feeding rate and a small exit rate the overall density is 
expected to be high. On the other hand for small a and large (3 the density is expected to 
be low. In addition, for a large system, where away from the boundaries the local density is 



expected to vary very slowly, the fundamental relation ([T0|) is expected to hold locally. Thus 



the current in the system cannot exceed a maximal current, as suggested by (|T0|). These 
features yield a rather rich phase diagram, as a and (5 are varied. 



We start by considering the phase diagram in the mean field approximation pT| , |22| , p^ 
Since correlations in this system are expected to be vanishingly small away from the bound- 
aries, this approximation is expected to yield a rather accurate phase diagram. In fact it 
turns out that the phase diagram obtained in this way is exact. 

To derive the mean field equations we note that the current Jj,j+i between sites i and 
2 + 1 is given by (o"t(l — 0".t+i)) for i = 1, . . . , — 1. In addition the currents at the two ends 
are given by Jq = — <^\) and Jat = (3{aN)- Neglecting correlations one finds that in the 
steady state, where all currents are equal, the following equations have to be satisfied: 

J = a{l- pi) = pi{l - P2) = . . . = pN-ii^ - Pn) = Ppn- (13) 

Solving these equations for J, pi, . . . , p^r the density profile in the steady state and the current 
are obtained. 

It is instructive to consider these equations in the continuum limit. Replacing pi by p{x) 
in ([TB| ) with < X < L yields the bulk current 

J(x)=pil-p)-D^, (14) 

where D is the diffusion constant which, by rescaling x, may be taken as 1. In this expression 
the first term represents the drive while the second term is the ordinary diffusion current. 
The evolution of the system is governed by the continuity equation dp/dt = —dJ/dx, 
together with the boundary conditions J(0) = a(l — p(0)) and J{L) = [3p{L). In the steady 
state the current J in (|1^ is independent of x yielding a density profile which has one of 
the two following forms 

p{x) = 0.5 + V tanh[t>(x — Xq)] 

p{x) = 0.5 + V coth[v{x — xq)], (15) 

where v"^ = 1/A — J. The two parameters xq and v (or alternatively the current J) are 
determined by the two boundary conditions, and are thus related to a and /?. 
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By matching the boundary conditions the density profiles and the current are obtained. 
The resulting phase diagram is given in Figure |^. The system is found to exhibit three 
distinct phases in the limit of large length L\ 

• Low density phase for which the bulk density is smaller than 0.5 with xq = 0{L). The 
density profile is basically fiat, except for a small region near the right end. In this 
phase p(0) = a and J = «(! — a). It exists for a < /? and a < 1/2. 

• High density phase with bulk density larger than 0.5 with xq = —0{L). The density 
profile is fiat except at a small region near the left end. Here p{L) = 1 — (3 and 
J = j3{l — (3). This phase exist for (3 < a and (3 < 1/2. The two phases co-exist on 
the line a = {3 < 1/2. 

• A maximal current phase in the region a > 1/2 and (3 > 1/2. Here the bulk den- 
sity is 1/2 exhibiting structures at both ends of the system. These structures decay 
algebraically as 1/a; when moving away from the ends. In this phase the current is 
maximal, namely J = 1/4. 

The phase diagram exhibits a first order line on which the high density and the low 
density phases coexist and two second order lines separating these phases from the maximal 
current phase. Typical schematic density profiles in the various phases are also given in 
Figure ||. 

In the mean field approximation fluctuations are neglected, and therefore by itself this 
analysis may not serve as a demonstration that phase transitions do take place in Id away 
from thermal equilibrium. In fact mean field approximation yields phase transitions in 
equilibrium Id systems, where they are known not to exist. It is therefore important to 
examine the role of fluctuations in this driven system and demonstrate that indeed the 
phase transitions found within the mean field approximation remain when fluctuations are 
taken into account. This has indeed been demonstrated for the TASEP A method 
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FIG. 2. The (a, f3) phase diagram of the totally asymmetric exclusion process exhibiting low 
density, high density and maximal current phases. The thick line represents a first order transition 
while thin lines correspond to continuous transitions. Schematic density profiles in the various 
phases are given. 

which goes beyond the mean field approximation and allows exact calculations of steady 
state properties of some driven Id models is described in the next section. 

IV. Matrix method 



A matrix method for calculating some steady state properties of the TASEP was intro- 



duced a few years ago The method has since then been generalised and applied to 

other models of driven systems. In the following we briefly outline the method as applied 
to the TASEP with open boundary conditions described above. 

We are interested in calculating the steady state distribution function P(cri, . . . , a^). In 
the matrix method one tries to express the distribution function by a matrix element of a 
product of particular matrices. Let D and E be two square matrices and {W\ and \V) be 
two vectors. For any given configuration, (ai, . . . , ctat) one considers the matrix product in 
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which each occupation number (Ji is replaced by either a matrix D or a matrix E depending 
on whether it is 1 or 0, respectively. The key question is whether one can find matrices 
D and E and vectors {W\ and \V) such that P(cri, . . . , a^) is proportional to the {W\, \V) 
matrix element of this product. Within this representation the distribution function may 
be written as 

N 

Piai, . . . , a^) cx (1^1 l[[a,D + (1 - a,)E]\V). (16) 

i=l 

A priori, it is not at all clear that such representation is available. However, if such represen- 
tation exists, it may yield a straightforward (though sometimes tedious) way for calculating 
the distribution function. For example, the density at, say, site i = 1 may be expressed as 



Pi = ^T.iW\DX2...Xr,\V), (17) 

{X,} 



where Xj = D, E, and the normalization factor is given by 

Zn=Y1 {W\XiX2 ...Xn\V) = (ly |C^|K), (18) 

{Xi} 

with 

C = D + E . (19) 

Thus, we may rewrite ( |I7D as 



pi = -^{WIDC'^-'IV) . (20) 

Densities at other sites and density- density correlation functions may similarly be expressed 
by other matrix elements. 

The main question at this point is how to find matrices and vectors such that (0) holds. 
To this end we consider the local currents in the system. Using the matrix representation, 
the current between sites i and i + 1 {i = 1, . . . , N — 1) may be expressed as 

Ji,i+i = -^{W\C'-'DEC''-'-'\V), (21) 
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and the currents at the two ends 



J, = ^{W\EC''-'\V) 
JM = ^mC''''D\V). (22) 

In the steady state all currents are equal. Taking matrices E and vectors {W\, \V) which 
satisfy 

DE = D + E (=C) 

a{W\E = {W\ (23) 
f3D\V) = \V), 

guarantee that all currents are equal, with 

J=^. (24) 



The question is whether these relations (p3D are sufficient to guarantee that the resulting 
distribution ([T6| ) is a steady state. For (0) to be a steady state one has to make sure that 

dP{ai,...,aN)/dt = 0, (25) 

for each of the 2^ microscopic configurations. The relations (|23|) only directly guarantee 
that + 1 of these 2^ equations are satisfied. This may suggest that (|23|) may not be 
sufficient to guarantee that the steady state distribution is given by (p!6D. However it can 
be shown, by direct inspection of Equations (pSj) that (p3D yields the steady state of the 
system. 

The problem is thus reduced to ffist finding matrices and vectors which satisfy (|16|), 
and then calculating some matrix elements to obtain, for example, the current J. It is 
straightforward to show that for a + (3 ^ 1 the matrices which satisfy (p!6| ) have to be of 
infinite order. Such matrices have been found and the current and density profiles and other 
correlation functions have been calculated E^ . The resulting phase diagram coincides with 



that obtained by the mean field approximation, although the density proffies are different. 
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For example, in the algebraic, maximal current, phase the local density decays to the bulk 
density like at large distances from the boundary, unlike the mean field result which 

yields a l/x profile. 

The matrix method proved to be very powerful in yielding steady state properties of 
TASEP dynamics. It has been applied and generalised to study partially asymmetric exclu- 
sion processes (ASEP) |]3l|-0^ and models with more than one type of particles In 
addition, replacing matrices by tensors proved to be useful in some cases [^,^]. However 
the method is restricted to one dimension. It is not standard in the sense that it can not be 
applied to any dynamical model. Moreover, there is no simple way to tell a priori whether 
or not it may be applicable for a specific model. 



V. Spontaneous symmetry breaking in one dimension 

In this section we consider a simple dynamical model which exhibits spontaneous sym- 
metry breaking (SSB) in Id. 

The model may be pictorially described in the following way: consider a narrow bridge 
connecting two roads. Cars travelling on the bridge in opposite directions do not block each 
other, although they may slow the traffic fiow in both directions. We assume that the two 
roads leading to the bridge from both sides are statistically identical. Namely the arrival 
rates of cars at the two ends of the bridge are the same. This system clearly has a right-left 
symmetry. Thus if this symmetry is not spontaneously broken one would expect that the 
long time average of the current of cars travelling to the right would be the same as that 
of cars travelling to the left. The question is whether the bridge is capable of exhibiting 
breaking of the right-left symmetry, and spontaneously turning itself into a 'one-way' street, 
where the current in one direction is larger than the current in the other direction. It turns 
out that this may indeed take place in the limit of a long bridge, demonstrating that SSB 
may take place in Id nonequilibrium systems. 

To model the 'bridge' problem we generalise the TASEP discussed in the previous section 
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37| , |40| ; ^ . We consider a Id lattice of length A^. Each lattice point may be occupied by 



either a (+) particle (positive charge) moving to the right, a (— ) particle (negative charge) 
moving to the left or by a vacancy (0). In addition positive (negative) charges are supplied 
at the left (right) end and are removed at the right (left) end of the system. 

The dynamics of the model is defined as follows: at each time step a pair of nearest 
neighbour sites is chosen and an exchange process is carried out 

+ 0^0+, 0-^-0, +-^-+, (26) 

with rates 1, 1 and q, respectively. Furthermore, at the two ends particles may be introduced 
or removed. At the left boundary (z = 1) the processes 

0^+ , -^0, (27) 

take place with rates a and (3, respectively. Similarly, at the right boundary (z = A^), one 
has the processes 

0^- , +^0, (28) 

with rates a and (3, respectively (see Figure ^. In the 'bridge' language the boundary terms 
may be viewed as traffic lights which control the feeding and exit rates at the two ends. 

Since the parameters a and f3 are the same on both ends of the systems the dynamics 
obviously possesses a right-left symmetry. The question of interest is whether or not this 
symmetry is preserved in the steady state. Clearly, for small a and large f3 the density of 
particles in the system is expected to be low, the two types of particles do not block each 
other, and the steady state is expected to be symmetric. On the other hand for (3 much 
smaller than a, particles are blocked in the system, the density is high and it is possible 
that symmetry breaking takes place. 

We start by considering the mean field approximation. It is straightforward to derive 
the mean field equations for the steady state. They take the form 

J+ =Pi[l- Pi+i - (1 - q)ini+i] 

J_ = mi+i[l - rrii - {1 - q)pi], (29) 
16 
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FIG. 3. A schematic representation of the input and output rates of the "bridge' model. 

for i = 1, . . . , N — 1, where pi and rrii are the densities of the (+) and (— ) particles at 
site i, respectively, and J+ and J_ are the currents of the positive and negative particles, 
respectively. In addition to the bulk equations one has four other equations for the currents 
at the boundaries 



J+ = a{l -pi- nil) = Ppn 
J_ = Pnii = a{l — Pn — ^n)- 



(30) 



These (2A^ + 2) equations may be solved numerically for (pi, . . . ,pjv; mi, . . . , m^v; J+, J-) to 
yield the {a, (3) phase diagram of the model. It is found that for large (3 the steady state 
is symmetric (with J+ = J J) while for small (3 the two currents are unequal in the steady 
state. 

The matrix method discussed in the previous section has been generalised and applied 
to this model [^|. However it turned out that a self-consistent matrix representation could 
be found for this model only for /5 = 1 or in the limit a oo. The limit a — oo is trivially 
mapped on the single species TASEP model. For /3 = 1 a phase transition is found although 
no spontaneous symmetry breaking takes place. According to mean field SSB is expected 
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FIG. 4. The time evolution of the current difference in the broken symmetry phase for 
{a = q = 1, f3 = 0.15, = 80). Time is measured in units of Monte Carlo sweeps. Flips be- 
tween the two symmetry related states are clealy seen. 

only at much lower exit rates (3. 

To demonstrate that the non-symmetric state found in the mean field approximation for 
small j3 survives fluctuations, numerical simulations of the dynamics have been carried out 
and the current difference J+ — J_ has been measured as a function of time for small (3, 
where the mean field approximation predicts a broken symmetry phase |^ . A typical time 
evolution of the current difference for a system of size = 80 is given in Figure ^. The 
figure suggests that the system flips between two macroscopic states: one with a positive net 
current and the other with negative net current. In the first case the system is predominantly 
loaded with positive charges moving to the right while in the second case it is loaded with 
negative charges moving to the left. This time course is characterised by a time-scale t{N) 
which measures the average time between flips. Clearly, when averaged over time, the 
current difference vanishes, yielding a symmetric state. This is to be expected since we are 
dealing with a finite system, and one certainly does not expect SSB to take place in a finite 
system. The question is how does the system behave in the thermodynamic limit N ^ oo, 
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000+ + + + + + 

X y 

FIG. 5. A typical microscopic configuration of the "bridge" in tlie broken symmetry state for 
small /?. Usually x is a small number of 0(1) and y is of 0{N) or vise versa. 

and particularly how does t{N) grow for large A^. Numerical simulations suggest that r 
grows exponentially with N. This means that the probability of a flip is negligibly small in 
a large system, and thus SSB takes place. 

In order to gain some insight into the flipping process we consider the limit of very small 
f3 pO| , ^ . In this limit, particles leave the system at a very small rate, and the system is 



filled with either positive charges moving to the right or negative charges moving to the left. 
Starting with a positively charged system, one would like to understand the mechanism by 
which a system of finite length flips into a negatively charged one. The evolution in the 
small f3 limit may be described as follows: with rate f3 a positive charge leaves the system at 
the right end. The vacancy created at this end may either move to the left with velocity 1 or 
may be filled with a negatively charged particle which in turn moves to the left with velocity 
q. When the negative charge reaches the other end of the system it is delayed for a while, 
but eventually leaves the system in time of order During this time, other negative 

charges may arrive at the left end forming a small blockage of negative charges. A typical 
configuration is given in Figure ^ It is composed of a segment of x negative charges at the 
left, another segment of y positive charges at the right and in between a segment of N—x—y 
vacancies. This configuration is denoted by {x,y). In the small /3 limit other configurations, 
for example those in which vacancies are present inside the charged segments, do not play a 
role in the global dynamics and may be neglected. 

The dynamics restricted to {x, y) configurations is rather simple. It may be viewed 
as the dynamics of a single particle diffusing on a square lattice performing the following 
elementary moves: 

{x,y) A (x+ 1,?/- 1) 
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{x,y) {x,y~ 1) 

{x,y)^{x-l,y+l) (31) 
{x,y) (x - l,y), 



where 

1 , a 



(32) 



2(1 + a) ' 2(1 + a)' 

are the rates of the various moves. Here the first two moves correspond to a positive charge 
leaving the system at the right end and being replaced by a negative charge (vacancy), 
respectively. Similarly the last two moves correspond to a negative charge leaving the 
system at the left end and replaced by a positive charge (vacancy), respectively. A schematic 
representation of this process is given in Figure |^. 

The biased diffusion process takes place as long as the particle stays within the triangle 
(x>0, y > , X + y < N). When it reaches the boundary of the triangle, for example 
X = the negative charge blockage at the left end disappears and on a very short time scale 
(as compared with 1/(3) the system is filled with positive charges from the left end, moving 
to (0, A^). The evolution of a positively charged system is thus represented by a random walk 



starting at (0, A^) with elementary steps defined by (pTj). Due to the bias of these elementary 
steps, a typical walk for large A^ ends on the x = axis. Once it reaches this axis it moves 
back to (0, A^) and the process starts again. This process repeats itself until the diffusing 
particle performs a walk which starts at (0, A^) and ends on the y = axis without touching 
the X = axis while diffusing. When this happens the blockage of positive charges at the 
right end is removed, the system is rapidly filled with negative charges moving to the other 
end of the triangle (A^, 0). This corresponds to a flip. The probability of such a walk taking 
place has been calculated, yielding the following flipping time [Q: 

r{N) = ^AT^/^e'^^, (33) 
P 

where C is a constant and 
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FIG. 6. A representation of the dynamics of the "bridge" in the hmit of smah /? as a biased 
diffusion process on a square lattice. The elementary moves and their rates are indicated. Starting 
with a positively charged system, namely at (A^, 0), possible trajectories are given. Typical walks, 
r , end on the x = axis. Walks, s , which end on the y = axis are rare, and they correspond to 
a flip. 
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K = In 



-(l-a - y^(l-a)(l-2a) )]. (34) 



The exponential flipping time is a direct result of the fact that the random walk correspond- 
ing to the evolution of the model is biased. It takes an exponentially long time to reach a 
distance of order against a bias. 

A very interesting question is related to the behaviour of nonequilibrium systems with 
spontaneous symmetry breaking when an external symmetry breaking field is introduced. 
In thermal equilibrium a symmetry breaking field makes the phase unfavoured by the field 
metastable or even unstable (when the field is large). For example, a positive magnetic field 
applied to an ordered ferromagnetic Ising system, removes the degeneracy between the two 
magnetic states. Only the state with positive net magnetization remains stable. The two 
magnetic states coexist only at zero field. This is a direct consequence of the Gibbs phase 
rule. 



It is known that in nonequilibrium systems, this is not necessarily the case p3| , p 
Namely when a symmetry breaking field is applied, the state unfavoured by the field may 
stay as a stable thermodynamic state. This is in violation of the Gibbs phase rule, which 
does not hold in nonequilibrium. 

The 'bridge' model described in this section provides a clear example for this behaviour 
40| . To demonstrate this point we introduce a symmetry breaking field by imposing bound- 



ary conditions which favour, say, the positively charged state, thus explicitly breaking the 
symmetry. More specifically, we consider an exclusion model where, instead of having bound- 
ary rates for both types of particles, we take for the positive charges and 
for the negative charges. The two exit rates are taken to be of the form (3± = /3(1 =F H), 
where < if < 1 is the symmetry breaking field, favouring the positively charged state. 
The analysis presented above in the limit /? — may be repeated for non-vanishing field H 
and the stability of the two phases may be analysed. Here again the dynamics is reduced 
to a diffusion process of the type (^) but with modifies rates (see Figure 0). Clearly the 
positively charged state is stable since it is favoured by the field. The question is whether 
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a(l-H) b(l-H) 

FIG. 7. The elementary moves of the biased diffusion process corresponding to the smah (5 
Hmit when a symmetry breaking field H is present. 

the negatively charged state is stable when the field H is non-vanishing. To examine this 
problem, we start with a negatively charged system (A^, 0) and consider a random walk de- 
fined by the diffusion process. It is easy to see that as long as if < a/(l — a) = 1/(1 + 2a) 
the walk is biased in the negative y direction, yielding a flipping time exponential in the 
system size. Thus the negatively charged state is stable even when it is unfavoured by the 
symmetry breaking field. 



VI. Pase separation in one dimension 

A phenomenon closely related to spontaneous symmetry breaking is that of phase sepa- 
ration. In Id equilibrium systems with short range interactions phase separation does not 
take place, and therefore no liquid-gas like transition is expected. The density of particles 
in such a system is thus macroscopically homogeneous. 

Recent studies have shown that driven systems may exhibit phase separation in Id even 
when the system is governed by local dynamics p9| , |45| -^. Several models have been intro- 



duced to demonstrate this behaviour. In these models more than one type of particles is 
needed for phase separation to take place. In the following we consider in some detail one of 
these models, and analyse the mechanism leading to nonequilibrium phase separation |39[|. 



The model is defined on a Id lattice of length with periodic boundary conditions. 
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Each site is occupied by either an A, B , or C particle. The evolution is governed by random 
sequential dynamics defined as follows: at each time step two neighbouring sites are chosen 
randomly and the particles of these sites are exchanged according to the following rates 



BC^ CB (35) 
CA^ AC. 



The rates are cychc in A, S and C and conserve the number of particles of each type A?^, Nb 
and Nc-i respectively. 

For q — 1 the particles undergo symmetric diffusion and the system is disordered. This is 
expected since this is an equilibrium steady state. However ior 1 the particle exchange 
rates are biased. We will show that in this case the system evolves into a phase separated 
state in the thermodynamic limit. 

To be specific we take g < 1, although the analysis may trivially be extended for any 
q ^ 1. In this case the bias drives, say, an A particle to move to the left inside a B 
domain, and to the right inside a C domain. Therefore, starting with an arbitrary initial 
configuration, the system reaches after a relatively short transient time a state of the type 
. . . AABBCCAAAB ... in which A, B and C domains are located to the right of C, A and 
B domains, respectively. Due to the bias g, the domain walls . . . AB . . ., ... BC . . ., and 
. . . CA . . ., are stable, and configurations of this type are long-lived. In fact, the domains 
in these configurations diffuse into each other and coarsen on a time scale of the order of 
q~\ where I is a typical domain size in the system. This leads to the growth of the typical 
domain size as (lni)/| lng|. Eventually the system phase separates into three domains of the 
different species of the form A. . . AB . . . BC . . .C . A finite system docs not stay in such a 
state indefinitely. For example, the A domain breaks up into smaller domains in a time of 
order q-'>^i'n{^B,NG} _ jj^ ^j^g thermodynamic limit, however, when the density of each type of 
particle is non vanishing, the time scale for the break up of extensive domains diverges and 
we expect the system to phase separate. Generically the system supports particle currents in 
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the steady state. This can be seen by considering, say, the A domain in the phase separated 
state. The rates at which an A particle traverses a B (C) domain to the right (left) is of the 
order of {q^^')- The net current is then of the order of —Q^^', vanishing exponentially 
with A^. This simple argument suggests that for the special case Na = Nb = Nc the current 
is zero for any system size. 

The special case of equal densities Na = Nb = Nc provide very interesting insight into 
the mechanism leading to phase separation. We thus consider it in some detail. Examining 
the dynamics for these densities, one finds that it obeys detailed balance with respect to some 
distribution function. Thus in this case the model is in fact in thermal equilibrium. It turns 
out however that although the dynamics of the model is local the effective Hamiltonian 
corresponding to the steady state distribution has long range interactions, and may thus 
lead to phase separation. This particular mechanism is specific for equal densities. However 
the dynamical argument for phase separation given above is more general, and is valid for 
unequal densities as well. 

In order to specify the distribution function for equal densities, we define a local occupa- 
tion variable {Xj} = {Ai, Bi,Ci}, where Ai, Bi and Cj are equal to one if site i is occupied 
by particle A, B oi C respectively and zero otherwise. The probability of finding the system 
in a configuration {Xj} is given by 

W^{{X,}) = (36) 

where 7i is the Hamiltonian 

n{{x,}) = E E (1 - -^)iCiB,+k + Ac+k + B,A+k) - (iv/3)^ (s?) 

i=l k=l 

and the partition sum is given by Zjy = q"^^^^^^^ ■ The value of the site index (z + k) in 



(p7|) is taken modulo N. In this Hamiltonian the interaction between particles is long range, 
growing linearly with the distance between the particles. 

In order to verify that the dynamics (B^) obeys detailed balance with respect to the 



distribution function ( PB| , |57D it is useful to note that the energy of a given configuration may 
be evaluated in an alternate way. Consider the fully phase separated state 
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A...AB ...BC ...C (38) 

The energy of this configuration is E = 0, and, together with its translationally relates 
configurations, they constitute the A^— fold degenerate ground state of the system. We now 
note that nearest neighbour (nn) exchanges AB BA, BC ^ CB and — > AC cost one 
unit of energy each, while the reverse exchanges result in an energy gain of one unit. The 
energy of an arbitrary configuration may thus be evaluated by starting with the ground state 
and performing nn exchanges until the configuration is reached, keeping track of the energy 
changes at each step of the way. This procedure for obtaining the energy is self consistent 
only when the densities of the three species are equal. To examine self consistency of this 
procedure consider, for example, the ground state (pH]), and move the leftmost particle A 
to the right by a series of nn exchanges until it reaches the right end of the system. Due 
to translational invariance, the resulting configuration should have the same energy as (pSD, 
namely E = 0. On the other hand the energy of the resulting configuration is E = Nb — Nc 
since any exchange with a B particle yields a cost of one unit while an exchange with a C 
particle yields a gain of one unit of energy. Therefore for self consistency the two densities 
Nb and Nc have to be equal, and similarly, they have to be equal to N^. 

The Hamiltonian (|37D may be used to calculate steady state averages corresponding to 
the dynamics (^). We start by an outline of the calculation of the free energy. Consider 
a ground state of the system (pSf ). The low lying excitations around this ground state 
are obtained by exchanging nn pairs of particles around each of the three domain walls. 
Let us first examine excitations which are localized around one of the walls, say, AB. An 
excitation can be formed by one or more B particles moving into the A domain (equivalently 
A particles moving into the B domain). A moving B particle may be considered as a walker. 
The energy of the system increases linearly with the distance traveled by the walker inside 
the A domain. An excitation of energy m at the AB boundary is formed by j walkers passing 
a total distance of m. Hence, the total number of states of energy m at the AB boundary 
is equal to the number of ways P{m) of partitioning an integer m into a sum of (positive) 



26 



integers. This and related functions have been extensively studied in the mathematical 
literature over many years. Although no explicit general formula for P{m) is available, its 
asymptotic form for large m is known 

P(m) ~ -^—^ (^(2/3)^^^ (39) 
Also, a well known result attributed to Euler yields the generating function 

oo 1 

Y=Y: q-'PM = -— , (40) 

m=0 {(Ijco 

where 

{q)^=lim{l-q){l-q')...{l-q-). (41) 

This result may be extended to obtain the partition sum Zn of the full model. In the limit of 
large N the three domain walls basically do not interact. It has been shown that excitations 



around the different domain boundaries contribute additively to the energy spectrum ^9 
As a result in the thermodynamic limit the partition sum takes the form 

Zn = iV/[(g)oo]^ (42) 

where the multiplicative factor N results from the A^— fold degeneracy of the ground state 
and the cubic power is related to the three independent excitation spectra associated with 
the three domain walls. 

It is of interest to note that the partition sum is linear and not exponential in A^, as is 
usually expected, meaning that the free energy is not extensive. This is a result of the long- 
range interaction in the Hamiltonian and the fact that the energy excitations are localized 
near the domain boundaries. 

Whether or not a system has long-range order in the steady state can be found by 
studying the decay of two-point density correlation functions. For example the probability 
of finding an A particle at site i and a B particle at site j is. 



{AB,} = ^ E (43) 
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where the summation is over all configurations {Xk} in which A'^ = A^^ = Nc- Due 
to symmetry many of the correlation functions will be the same, for example (AiAj) = 
(BiBj) = (CiCj). A sufficient condition for the existence of phase separation is 

lim Mm ((AiAr) ~ {Ai)(Ar)) > 0. (44) 

Since (Ai) = 1/3 we wish to show that limr^oo MmN^oo{AiAr) > 1/9. In fact it can be 
shown that for any given r and for sufficiently large A^, 

{A,Ar) = 1/3 -0{r/N). (45) 

This result not only demonstrates that there is phase separation, but also that each of the 
domains is pure. Namely the probability of finding a particle a large distance inside a domain 
of particles of another type is vanishingly small in the thermodynamic limit. 

Numerical simulations of the model for the case of unequal densities, where such analysis 
cannot be carried out, strongly indicate that phase separation takes place as long as none of 
the three densities vanish. They also indicate that the coarsening process which accompanies 
phase separation is rather slow, with the characteristic length diverging like In t at long times. 



VII. Summary 

In these lecture notes some collective phenomena which occur in one-dimensional driven 
systems have been reviewed. These systems have been extensively studied in recent years 
by introducing simple models and analysing their steady state properties. Some of these 
models have been demonstrated to exhibit a rich variety of phenomena which are unexpected 
in equilibrium one- dimensional systems. 

Simple asymmetric exclusion processes in open systems were shown to exhibit both 
first order and continuous phase transitions. Other systems which have in the past been 
demonstrated to exhibit phase transitions in Id are directed percolation |^9|,Q and contact 



processes |5T]]. These system, however, possess one or more absorbing states. Once the 
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system evolves into one of these states the dynamics is such that the system is unable to 
exit. Under these conditions, the existence of a phase transition between a trapped and an 
untrapped states is rather natural. Usually, once the dynamics in these models is generalised 
to allow for an exit from the absorbing state no phase transition takes place. The phase 
transitions occurring in the asymmetric exclusion processes discussed in this paper are rather 
different, as the dynamics in these models does not possess absorbing states. 

Mechanisms which lead to spontaneous symmetry breaking and phase separation in one- 
dimensional nonequilibrium systems have been discussed. A common crucial feature of these 
models is that the dynamics conserves, at least to some degree, the order parameter. In the 
'bridge' model, the densities of the two types of particles are conserved in the bulk although 
they are not conserved at the two ends of the system. In the ABC model, on the other hand, 
the three densities are fully conserved. When non-conserving processes are introduced into 
these models spontaneous symmetry breaking and phase separation do not take place. It 
would be very interesting to consider the possibility of spontaneous symmetry breaking in 
one dimension when the dynamics does not conserve the order parameter. A related problem 
has been considered in the context or error correcting computation algorithems. An example 
of a one- dimensional array of coupled probabilistic cellular automata has been constructed 
and shown to yield breaking of ergodicity, as would a model with spontaneous symmetry 
breaking This approach suggests that indeed spontaneous symmetry breaking in Id 

may exist even when the dynamics is not conserving. However the example given is rather 
complicated and not well understood. 

In spite of the progress made in recent years in the understanding of nonequilibrium 
collective phenomena, many basic questions remain open, even for the restricted and rela- 
tively simple class of systems which evolve into a steady state. For example a classification 
of continuous nonequilibrium transitions into universality classes, like the one which exists 
for equilibrium transitions, is not available. Also the djTiamical process of the approach to 
steady state is far from being understood in many cases. The approach outlined in these 
notes, which involves constructing simple dynamical models and analysing the resulting col- 
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lective behaviour may prove to be helpful in developing better understanding of some of 
these complex questions. 
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